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ANNUAL  PROGRESS  REPORT 

GRANT  #:  N00014-89-j-3001  R&T  Code:  4412068 

PRINCIPAL  INVESTIGATOR:  Sangtae  Kim 
INSTITUTE:  University  of  Wisconsin  -  Madison 

GRANT  TITLE:  Microstructural  Models  of  Interactions  that  Govern  Protein  Conforma¬ 
tions 

PERIOD  OF  PERFORMANCE:  1  July  1989  -  31  May  1990 

OBJECTIVE:  To  apply  computational  methods  for  microstructures  suspended  in  a  viscous 
fluid  to  the  calculation  of  the  slower  dynamic  processes  in  protein  folding.  The  concerted 
motion  of  the  entire  tertiary  structure  will  be  considered  with  the  secondary  consituents 
modeled  as  sohd  segments.  The  short  term  objective  is  to  determine  diffusion/ffiction 
coeflicients,  first,  for  the  overall  protein  and  then  for  the  larger  sub-structures. 

ACCOMPLISHMENTS  (last  12  months):  My  students  and  I  have  completed  the  inter¬ 
facing  of  the  Brookhaven  Protein  Database  to  our  simulation  code  for  the  hydrodynamic 
mobility/friction  factor,  the  Completed  Double  Layer  Boundary  Integral  Equation  Method 
(CDL-BIEM).  The  key  step  wais  accomplished  by  converting  M.L.  Connolly’s  Moleculair 
Suiface  Package  (Piecewise  Quartic  Molecular  Surface,  Surfa.ee  Rendering  by  Foliation, 
Triangulation  by  Recursive  Bisection)  into  a  boundary  element  mesh  generator  for  CDL- 
BIEM.  This  process  creates  realistic  three-dimensionaJ  rendering  of  surfaces  of  smaller 
protein  residues. using  about  300  boundary  elements.  An  entire  enzyme  such  as  lysozyme 
(Figure  1)  cam  be  rendered  using  about  3000  boundauy  elements. 

For  N  boundary  elements,  the  CDL-BIEM  calculation  involves  the  iterative  solution 
of  a  dN  by  3iV  dense  linear  system,  as  discussed  in  o\ir  publication  [1].  For  compact, 
sphere-like  shapes,  the  Jacobi  iteration  of  the  form,  x  =  Mx  -f-  6,  converge  very  quickly, 
in  as  few  as  five  iterations.  For  more  complex  shapes  associated  with  tertiary  structures 
of  globular  proteins,  the  procedure  converges  too  slowly.  We  have  implemented  more 
powerful  iterative  methods  based  on  the  residual  reduction  methods  of  Eisenstat  ti  al. 
(SIAM  J.  Numerical  Analysis  20,  pp.  231-154,  1983).  Typical  convergence  behavior  of  the 
mobility  coeflicient  (which  is  directly  proportional  to  the  diffusion  coefficient)  is  shown  in 
Figure  2  for  a  fairly  difficult  but  representative  problem  from  the  viewpoint  of  numerical 
analysis:  a  composite  object  formed  by  joining  a  thin  rod  with  a  sphere. 

The  largest  system  attempted  to  date  involved  the  calculation  of  the  mobility /diffusion 
coefficient  (or  equivalently,  the  friction  factor)  of  lysozyme,  using  the  shape  obtained  from 
the  Brookhaven  protein  database.  With  about  3000  triangular  boundsury  elements,  the 
resulting  linear  system  of  equations  is  9000  by  9000.  Using  a  memory-saving  version  in 
which  the  matrix  elements  are  recalculated  as  needed,  an  iteration  tsikes  about  10  hours  on 
a  super-minicomputer  nmning  at  10  MFlops,  and  a  complete  calculation  requires  about 
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10  iterations.  In  other  words,  a  Gigaflop  computer  will  be  able  to  perform  the  entire 
calculation  in  about  1  hour. 

SIGNIFIC AN CE;  It  is  well  known  that  there  is  a  systematic  pattern  of  deviation  between 
the  theoretical  predictions  and  experimental  measurements  of  diffusion  coefficients  as  a 
function  of  protein  shape.  The  more  simple,  spherical  shapes  are  adequately  explained 
by  theories  based  on  the  hydrodynamic  radius  and  the  radius  of  gyration.  The  ‘outliers’ 
are  invariably  those  proteins  with  contorted  surface  shapes.  The  CDL-BIEM  calculations 
show  great  promise  in  predicting  both  translational  and  rotational  diffusion  coefficients  of 
these  difficult  proteins  from  first  principles. 

The  significance  of  predicting  protein  diffusion  coefficients  using  CDL-BIEM  lies  in  the 
fact  that  we  will  be  able  to  test  the  Stokes-Einstein  theory  without  making  any  assumptions 
on  how  to  calculate  the  mobihty  of  the  rigid  bodies  we  use  to  model  the  proteins.  Past 
crdculations  of  diffusion  coefficients  have  reUed  on  three  sets  of  assumptions: 

1.  Assumptions  associated  with  the  generalized  Stokes-Einstein  theory. 

2.  Assumptions  on  how  to  generate  an  appropriate  molecular  surface  (using  the  given 
atomic  coordinates)  to  use  in  the  theory. 

3.  Assiunptions  on  how  to  calculate  the  mobility  of  the  rigid  body  defined  by  the  molec¬ 
ular  surface. 

We  can  eliminate  the  third  set  of  assumptions,  thus  shedding  light  on  the  first  two 
sets.  Since  the  first  two  are  also  involved  in  the  more  complicated  problem  of  protein 
folding  (or  dynamics  of  internal  motions)  we  gain  information  of  use  in  future  calculation 
of  these  motions. 


WORK  PLAN  (next  12  months):  The  objective  next  year  is  to  move  our  work  to  a  new 
generation,  GigaFlop  computer  architecture,  featuring  many  processors  (supercomputers- 
on-a-chip)  running  in  parallel,  to  be  located  in  the  new  Center  for  Parallel  Computationad 
Engineering,  UW-Madison.  Our  iterative  methods  are  highly  parallel  and  should  achieve 
close  to  the  theoretical  speedups  in  a  multi-processor  environment.  Runs  on  the  IBM  3090, 
Sequent  Symmetry  and  the  Intel  iPSC/860  have  all  confirmed  this  expectation. 

PUBLICATIONS  TO  DATE: 


1.  A  theoretical  paper  (submitted  before  the  contract  was  awarded)  heis  appeaired. 

Karrila,  S.J.,  Fuentes,  Y.O.  amd  Kim,  S.  (1989)  Parallel  computationad  strategies 
for  hydrodynamic  interactions  between  rigid  paurticles  of  aurbitrary  shape  in  a 
viscous  fluid,  J.  of  Rheology  33,  913-947  (copy  enclosed). 

2.  An  abstract  was  recently  submitted  to  the  session  on  Molecular  Biophysics  of  Proteins, 
Peptides  and  Polynucleotides,  chaired  by  M.L.  Yaumush,  in  the  1990  national  meeting 
of  the  Americam  Institute  of  Chemical  Engineers  (copy  enclosed). 

Kim,  S.  amd  Brime,  D.A.  (1990)  Theoretical  predictions  of  tramslationad  amd  ro¬ 
tational  diffusion  coefficients  of  proteins.  Abstract. 


□  □ 


i 


Figure  1.  Triangulation  of  lysozyme  for  the  boundary  element  calculation  of  the  friction 
coefficients. 


Figure  2.  Convergence  of  the  generalized  conjugate  residual  method  of  Eisenstat  et  al.  as 
compared  to  the  conjugate  gradient  method  applied  to  the  normal  equations. 
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